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Abstract The covariant generalizations of the background dark sector coupling suggested 
in G. Mangano, G. Miele and V. Pettorino, Mod. Phys. Lett. A 18, 831 (2003) are considered. 
The evolution of perturbations is studied with detailed attention to interaction rate that is 
proportional to the product of dark matter and dark energy densities. It is shown that some 
classes of models with coupling of this type do not suffer from early time instabilities in 
strong coupling regime. 
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1 Introduction 

Recent observations have shown that the Universe is spatially flat and presently accelerat- 
ing. Most attempts to explain this acceleration involve the introduction of dark energy, as 
a source in the Einstein field equations. In addition, the important contribution to the total 
density is the dark matter, whose existence is inferred indirectly by observing its gravita- 
tional influence on normal matter, such as stars, gas and dust. Standard ACDM model, in 
which dark energy is considered as a small positive cosmological constant and dark matter 
treated as gas of cold non-baryonic particles, provides a very good fit to the supernovae data 
[ 1 1 as well as CMB measurements [2.1 and observations of large scale structure (3) , but the 
small and fine-tuned value of the cosmological constant cannot be explained within current 
particle physics |4j. As a result a lot of other more expected from theoretical point of view 
cosmological models have been proposed to giving a dynamical origin to dark energy. 

Models with non-minimal coupling of dark matter and dark energy have called attention 
in the last decade. The nature of both dark energy and dark matter are still unknown, and 
in principle, the additional interaction between them is possible. Moreover, the coupling be- 
tween matter and quintessence is motivated by high energy particle physics considerations 
(see (5) for a set of references). There are a large number of works in where the inter- 
action is seen in models with quintessence 15 6 7 8 91011121 field. Also available the 
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phenomenological approach, in which the interaction is introduced into the conservation 
equations for dark matter (c) and dark energy (x), considered as perfect fluids. The back- 
ground energy exchange in the dark sector can be represented by 

p' c = -3,^p c + aQ, (1) 

p' x = -3Jff{l+w x )p x -aQ, (2) 

where the prime denotes derivative with respect to conformal time X, the bars mark the 
background quantities, a is the scale factor, Jif = a! /a, w x = P x /p x and the coupling Q is 
some function of background variables. 

In the simplest models the quantity Q is a linear combination of the dark sector densities 

Q=A c p c +A x p X7 (3) 

where Aj = 3a;H or A/ = J}. Here «/ are the dimensionless constants and i] are constant 
interaction rates 1131 . There are a large number of works in which for coupling <[3j) was 
considered the background [ 13 , 14, 15] and the perturbations evolution 1 16 , 17 , 18 19 ,20,21]. 

In the paper |22| in the context of the stimulated decay of dark energy into dark matter 
had been proposed the interaction of the form 

Q = YPc~pL (4) 

where a, /3 are the dimensionless constants, and y is a constant parameter with dimension 
[density '^"^ xtime -1 ]. This coupling provide a mechanism to alleviate the cosmic coin- 
cidence problem 1231 . [24|, namely why the dark energy density nearly coincides with the 
dark matter density presently. The models with interaction l[4} also have been studied in the 
framework of holographic dark energy [25]. By assuming that a = l,jS=0ora = 0,jS = l 
one obtains the particular cases of coupling Q. 

Although the basic properties of background solutions for the interaction (O are previ- 
ously considered, the perturbation evolution in such scenarios is not examined except the 
mentioned above special cases. In this paper we shall cover the perturbation equations and 
investigate the solutions of them with detailed attention to the most physically reasonable 
choice of parameters a = J3 = 1. 



2 Background. 

We consider the spatially flat Universe filled with radiation (r), massless neutrino (v), bary- 
onic matter (b), cold dark matter (c) and dark energy (x). We choose the linear parametriza- 
tion of the dark energy equation of state 1261 

w x {a) = w + wi(l -a) (5) 

with constants wq> — I and < w\ < —wo, so the dark energy is subdominant with respect 
to the radiation at early times and the phantom divide line w x = — I crossing does not occur. 
Since w' s /Jt p = — awi, in the radiation dominated era the dark energy equation of state 
parameter w x is well approximated by a constant wq + w\. 

Background dynamics in the presence of coupling I0 is completely described by Fried- 
mann's equation 

^2 = 8*G fl2 . (6) 
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Fig. 1 Background energy densities in units of po cr for a = ft = 1, Ho = 70 km s -1 Mpc _1 , £2 X = 0.70, 
£2 C = 0.2538, = 0.0462. Examples of dark sector densities evolution at negative and positive coupling are 
shown on left and right panels, respectively. 



and continuity equations 

p' r = -4Jfp r , p' v = -4J%>p v , (7) 

p' b = -3Jfp b , (8) 

p' c = -3,^p c -ayp?p?, (9) 

p' x = -3 JT( 1 + w x )p x + ayp?pP . (10) 

To determines the characteristic strength of the dark sector coupling it is convenient to 
introduce the dimensionless parameter A defined by 1 25 ] 

a = ip£*-V (id 

where po cr = 3Hq/(87vG) is the current critical density. 

System of equations (f6t-<l 10b can be numerically solved, as was done in [22] for a = 1 
and the dark energy equation of state w x = — 1. However, the equations I0 and d 10b imme- 
diately suggest several analytical conclusions. 

When y < 0, the dark energy density drops with the scale factor «(t) not slower than 
a -3(i+n ) an( j the dark matter density can increase under laypc ^Px /«^H > 3 condition. 
These features offer a simple way out to the cosmic coincidence problem, since in this 
models the present situation of p x ~ p c have occurred many times in the past 1 22 1 . 

The class of models with positive y also describe a interesting scenarios, as it allow to 
consider the dark energy grew out due the dark matter decay. However, the interaction of this 
type potentially leads to unphysical results, because at early times the dark energy density 
could be driven to negative value in some models with decaying dark matter 1 16 1. 

Figure 1 shows the time evolution of the dark matter and dark energy densities at values 
parameters a = jS = 1 . The pattern of the solutions varies depending on the model parame- 
ters. When w x < — 1 /3 at high redshifts and y < 0, one can obtain a sequence of regimes of 
the weak 

max(\ayp c /J^\,\ayp x /Jf\) < 1 and strong max(\ayp c /Jtf'\ , \ayp x /Jf |) > 1 dark sector 
coupling. At the stage of the weak interaction under the condition w x < — 1 /3 the dark en- 
ergy density decays more slowly than a~ 2 (f), and hence in the radiation dominated era the 
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inequality \ayp x /Jif\ > 3 will be reached at some time, from which the density of dark mat- 
ter increases. The increased dark matter density increases the dark energy decay rate, the 
condition \ayp x /Ji?\ <C 1 recovers and dark energy has no longer significant impact on the 
evolution of dark matter. Moreover, since the dark matter density is now reduced propor- 
tional to the inverse third degree of the scale factor, there comes a point when the influence 
of dark matter to dark energy also becomes negligible, i.e. a next stage of a weak interaction 
begins. For — 1/3 < w x however, if the interaction is weak initially, it will remain so through 
the subsequent evolution. This means that the stage of strong interaction in such models is 
only one. Accordingly, in such models the density of dark matter in the radiation dominated 
era is extremely small. 

The right panel of Figure [TJ depicts the time evolution of the density of dark matter and 
dark energy at values parameters a = jS = 1 and positive y. In this case, the density of dark 
energy is always positive. 



3 Model and perturbed equations. 

Let us consider a spatially flat FLRW Universe with small scalar type perturbations. The 
line element may be written as 

ds 2 = a 2 (T) {-(1 +2(j))dT 2 +2BjdTdx i + [(1 - 2\j/)8ij + 2/:.,, >/av/v ; J . (12) 

The energy-momentum tensor of a perfect A-fluid is given by 

T£ v = (p A + P A )i£ HA y+P A 8?, (13) 

where p A = Pa(1 + 8 A ) is density, P A = P A + 8P A is pressure, and the fluid four- velocity 
= dx^/ds at linear order of the perturbations is 

«£ = ^[(W),vi], «v=fl[-(l + 0), VA.i + Bj}. (14) 

We assume that the anisotropic stress of dark energy and dark matter vanishes, so the 
dark species can be treated as perfect fluids. In the general case of coupling fluids divergence 
of the energy-momentum tensor of each fluid yields 

where 4-vectors are restricted by the constraint 1271 

Ie£=o, (i6) 

A 

which follow from the conservation law of the total energy-momentum tensor. For conve- 
nience one can decompose these 4-vectors into two parts 

QH = QA«! i +F£, Qa = Qa + 8Q a , M M Ff = 0. (17) 

Here is the total four- velocity, Q A is the energy density transfer rate and is the momen- 
tum density transfer rate of A-fluid in the total matter gauge. By definition, =aT x (0, f A ), 
where f A is a momentum transfer potential 1161 . 
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The energy exchange in the background does not determine the covariant form of energy 
exchange. Instead, an energy exchange four-vectors Q% must be specified. The simplest 
scalars, which can be used to construct quantities Qc , Q x are 

T c %=p c , T x %=p x -3p x , T c %T x v a =p c p x (18) 

where the last relation holds at linear order. Accordingly, the direct generalization I0 with 

G t . = -& = -yp?p!!, (i9) 

is covariant at the first order. 

We are interested first of all in models with decaying dark matter or with decaying dark 
energy. The most natural way is to choose a vector F A in the form 

F» = -F» = ~ (0, (ypfpf (v - v c )Y) (20) 

or 

= -Ft 1 = ~ (0, (YP?p! (v - v x )y) . (21) 

The momentum transfer vanishes in the dark matter rest frame or in the dark energy rest 
frame respectively. A few more general form of coupling can be written as 

G.v = -Q? = P?p!M + (22) 

where y c and y x are arbitrary constants. 

The continuity equations for coupled perfect fluids can be obtained by linearization of 
equations H51 . As a result, one can write 1161 

8p' A + 3Jt {8p A + 8P A ) - 3 {Pa+Pa) V - k 2 (Pa +Pa) (v A +E') = aQ A <i>+a8Q A , (23) 

[(Pa + P A )(va+B)}' + 4Jf(p A + P A ) (v A +B) + (p A + P A )$ + 8P A = aQ A (v + B)+ af A , 

(24) 

where for coupling (122b 

8Qc = SQ X = -(Yc+Yx)P?p!!(aS c +p8 x ), (25) 

fc = -fx = Pep! (Yc(y - v c ) + y x (v - v x )) . (26) 

Pressure perturbations in the general case can be expressed in terms of the density per- 
turbations and the velocity potentials 

SP A = c] A 8 P A + (c 2 sA - c] A(ad) )(3,X?(l+ WA )p A - a Q A ) ^ . (27) 

where 9 A = —k 2 (B + v A ), c^/^ = P A /p A is the adiabatic sound speed and the rest frame 
sound speed s s a defined by 



8P A 
8 Pa 



(28) 

A— fluid rest frame 



We will work in synchronous orthogonal gauge <j> = 0, B = in the Fourier space, using 
the notation of 1281 : 

% = TJ, k 1 E k = -h/2-3r 1 . (29) 
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Now the conservation equations 015b for the dark energy and dark matter in the syn- 
chronous gauge takes the form 

8' x + 3Je(c 2 sx - w x )8 x + (1 + w x ) d x + 3 J? [3Jf(l + w x ) (4 - w x ) + w' x ] j r 



k 2 



1 +w Xjl 



a8 c + (P-l)8 x + 3J?(ci-w x )^ 



ti = { 7c + 7x )ap?pl- 1 

S' c + l -h' + 6 c = (Y c + Y x )apr'p! [(\-a)8 c -R8 x \, 



,1 2 



aril 8 - 1 



(30) 
(31) 

(32) 
(33) 

Note that in this phenomenological approach the sound speed c sx is needed to be fixed by 
hand [291. In the case of quintessence dark energy one have to set c sx = 1, and we adopt this 
value in the following. 

The perturbed Einstein equations are well known, and can be found in |28|. We repro- 
duce here only one of them 



B > + ^ (i - 34) e x - ff-s x = (Yc(e c - e x ) - (y c + y x y sx e x ) , 

1 + W x 1 -t- w x 

d^ + Jf0 c = Y x ap?- 1 pg[0 c -d x \. 



h" + Jfti = -8rcGV (5p + 38P) 



(34) 



Equations OOt - d33 b provide a set of coupled equations covering the dark sector density 
evolution. 



4 Large scale perturbations. 

The coupling terms appearing in the dark energy pressure perturbations may lead to the early 
time instabilities, as was first pointed out in [16]. The similar phenomenon is well known 
in inflationary multi-fields models, where on large scales entropy perturbations can source 
adiabatic ones [30]. To study this phenomenon here we write the second order differential 
equations for the dark energy density perturbations. This approach [18] allows to identify 
areas of possible instabilities before solving the perturbation equations. 

The perturbed fluid equations ( 130b - I l32b and equation d34b can be combined to 

S' x ' + A X ,?P8' X + B x ,3f z 8 X = C x Ji? 2 , (35) 

where on large scales (k <C Jf?) 

= 1 - 3wt - + {tt^: - (/3 - 1)r+ ) — /j (36) 

M" reap?/*!!- 1 \ 3w x w[. 
B x = 3{\-w x ) l--T=~+- 



jt? 2 (i+w x )je J {\+w x )je 

+(B n r am *~ l (2 1 rc+r + *P?p!- 1 . f 
—y 2+ j?2- 1+Wx +f 
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C r --(l+w,) 3 + -^ — __+/■_-_ 



3Jt?\3-3w x + 



_„ _8-l \ -(v-6-1 ,, 



{l+w x )je l+w x 



7 / «A : pr ^ / , ^' Yc + Y+aP?Px \ A y+ a PcPx ' a8 



aS c - 2 + 



J? 2 l+w x Jt? 



k 2 



Yja 2 pr- l P ?- 1 

JV 1 



+(l-a) 



8P 



-aS c +-{l+w x )[8 + 3— ), (38) 



and 



ln|«l -»|3 + 3w,-^^) + 5 



. 7+ = 7f + 7v (39) 



In case of minimal coupling and a constant dark energy equation of state the quantity / 
is zero at all times. When y f = or y A = 0, terms with 8 C can be ignored, since in the latter 
case one can work in particular synchronous orthogonal gauge in which the dark matter fluid 
has a vanishing velocity. 

The equation Q5b is written in slightly different form then was discussed in the Ref. [ 18 1. 
In particular, it involves h' instead of the time derivative of dark matter density perturbations. 
This difference can be valuable, since both quantities h', 8' c are related by equation Bil l and, 
in the general case, may implicitly depend on the dark energy perturbations. Used here 
notation is convenient to study of generations of non-adiabatic perturbations in the radiation 
dominated era. 

At first approximation, the source terms can be calculated using usual adiabatic mode 
solutions. Since the dark species are subdominant in early Universe, the corresponding con- 
tributions to h' and total 8, 8P are negligible. Actually, this is the assumption about the 
initial conditions, but it is certainly satisfied if the initial adiabatic conditions are imposed 
at the stage of weak coupling. This assumption breaks down if the dark energy perturba- 
tions have increased dramatically. It means that the right hand side of the equation d35 b can 
be treated at early times and on initial stages of dark energy inhomogeneities growth as an 
external force that is independent of the dark energy perturbations. In this approach, the 
negative sign of A x , B x or both of them indicate on the existence of large scale instabilities 
due a anti-damping force or exponential growth of intrinsic dark energy perturbations. The 
nature of instability can be revealed by considering the time evolution of the gauge-invariant 
curvature perturbation on uniform density hypersurfaces 



P 



(40) 



that is conserved on large scales for adiabatic perturbations. 

We carry out the detailed analysis of perturbations for two cases: y c 
Yc = 0, Yx = 7 with a = jS = 1, i.e. for 



y, Yx = and 



(41) 



where b = 1 for the first case and b = for the second. 
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For coupling above 

B. - 3(1 - ,,» (■ - £ + "/) - ^ + (43, 

Using equations ©, (9), d 1 Ob and d39b . the numerical values of these coefficients can be 
calculated directly from the background solutions. 

When y > 0, in the radiation dominated era one can assume <C 1, 3> 1 (see 
the right panel of Figure [T}. In this limiting case the dark energy has no effect on the evolu- 
tion of other fractions and their perturbations. By applying equation quantity / can be 
approximated as 

/«-2--^---^. (44) 
Keeping only the dominant terms, the equation d35b reduced to 



,b + 1 yap c s , , , ^ 2 _ , ( 2> yap c \ s _ ^ 2 2> + 1 / yap, 



2 



Coefficients A x and 5 V are positive and catastrophic growth perturbations do not occurs. In 
order of magnitude, equation ( 145b implies the estimation 

(46) 

For example, the "standard" adiabatic condition Sp c /p' c = 8p x /p[ taking into account the 
background equations ll9ll.lll0t in the radiation dominated era yields 

S, = -ffs c ~T (47) 

as for the adiabatic mode in the synchronous gauge 8 C T 2 . 

When 7 < 0, the coefficients A x , B x can take large negative values in regime of very 
strong coupling if \yap c / M'\ 3> 1. In scenarios with b = 1 and w x (a) < —1/3 at early times 
they are becoming both negative together under this condition . Hence such models suffer 
from the fast growth of non-adiabatic perturbations in early Universe. However, constraints 
on coupling with b = can be weakened, because in this case in some range of parameters 
wo, Wi the coefficient B x remains positive and increases in the strong coupling regime. 

To verify the analytical conclusions we have modified the public available CAMB code 
1311 . The initial adiabatic conditions for all non-dark species are imposed the same as in the 
non-interacting case. The initial values of dark sector perturbations are taken in accordance 
with 8p c /p' c = 8p x / p' x = 8p r /p' r . In the presence y x ^ the cold dark mater rest frame and 
synchronous frame are not coincide and it is not possible to adopt the CAMB conventions 
d c = consistently. In our numerical calculations the residual synchronous gauge freedom 
was fixed by choosing d c (Zi„) = 0, where T,„ is the time moment of the initial conditions 
setting. The results are shown in Figure[2] 
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a. = -ixio" 8 a,=-ixio" 




a 



Fig. 2 a) The curvature perturbation £ on super-Hubble scale in models with coupling j4U and constant dark 
energy equation of state w x < —1/3, Growth of the non-adiabatic perturbations begins at the stage of the 
strong interaction. b)The density perturbation evolution in model with w x > — 1/3 in radiation dominated era. 
Initial dark sector conditions are set at T,„ = 2 x 10~ 8 Mpc by Sp c /p' c = Sp x / p' x = Sp,-/p' n d x (?in) = 0- In 
both cases the background densities corresponds to present values Q x — 0.70, £2 C = 0.2538, Qi, = 0.0462, the 
initial curvature perturbation £ (t,„) equal to 1 x 10~ 25 and comoving wave number is k = 1 x 10~ 5 Mpc~' . 



5 Special case. 

The special case of I l22t is the interaction 

Q M =YPcp x {"c-^)- (48) 

In the linear perturbation theory the velocity potential is a first order quantity, and hence the 
spatial components of 4- vector Q 11 are proportional to the product of the densities of both 
dark species and relative velocity. In components 

G° = 0, Q^^ivi-vi). (49) 
a 

The background dynamics here is the same as at minimal coupling, but the perturbation 
evolution is different. The coefficients of second-order equation J35b for the large scale 
perturbations are 

2M" yap c 

Ax = l- 3Wx - — + J —^— -f, (50) 

j lx = 3^(l_ w ,)^l___ / j___ + __j. (51) 

Since now p c = p c oa~ 3 , we have \yap c jjff\ ~> 1 in early Universe. Thus, at y < the coef- 
ficients A x , B x are both negative and there is a rapid growth of non-adiabatic perturbations 
in the radiation dominated era. At positive y the non-adiabatic growth of long-wavelength 
perturbations do not occurs. These analytical results are confirmed by numerical computa- 
tions. 

In the short-wave approximation the equations d30b - d34l > can be combined into 
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^ 2 (K 5+3 y) + (OT( 3(1 -^ + J))- (53) 

At positive y all coefficients in the left hand side of these equations are also positive, 
what excludes the presence of small scale adiabatic instabilities. For instance, at early times 
in the strong coupling regime with ^fjif 2 ^> k 2 3> J^ 2 the first equation gives at leading 
order 

S' x = {l+w x )S' c . (54) 

Equations d52b and J53b take the same form as in the minimal coupling case under condi- 
tions c -^- <C 1 and c ^p- <C 1 respectively. In particular, when y — > 0, the second one reduced 
to the standard growth equation. 

6 Conclusions. 

We examined the covariant generalization of the coupling ©. The evolution of perturbations 
is studied paying particular attention to the most favored interaction rate that is proportional 
to the product of dark matter and dark energy densities. It is shown that the models of 
the form fl41b with b = 1 and w a (t) < —1/3 in radiation dominated era suffers from early 
time instabilities due fast growth of large scale non-adiabatic perturbations. Models with 
w x (t) > —1/3 in radiation dominated Universe are free from this defect. Also an interesting 
coupling d48b with positive y is viable. 
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